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Abstract 



In this paper the averaged weak (AWEC) and averaged null (ANEC) energy 
conditions, together with uncertainty principle-type restrictions on negative en- 
ergy ("quantum inequalities"), are examined in the context of evaporating black 
hole backgrounds in both two and four dimensions. In particular, integrals over 
only half-geodesics are studied. We determine the regions of the spacetime in 
which the averaged energy conditions are violated. In all cases where these 
conditions fail, there appear to be quantum inequalities which bound the mag- 
nitude and extent of the negative energy, and hence the degree of the violation. 
The possible relevance of these results for the validity of singularity theorems 
in evaporating black hole spacetimes is discussed. 
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1 Introduction 



It is by now well-known that quantum field theory permits violations of all of the local 
energy conditions used in classical general relativity These conditions are employed 
in a variety of ways, such as in the proofs of theorems on the occurrence of singularities 
in gravitational collapse and cosmology or of cosmic censorship. Two of the weakest 
such local conditions are the "weak energy condition" and the "null energy condition" 
which state that: 

T^^^>0, (1) 

for all timelike vectors u^, and 

T^K^K^O, (2) 

for all null vectors K^, respectively JJJ. Note that the null energy condition follows by 
continuity if the weak energy condition holds. These two local conditions are satisfied 
by known forms of classical matter, but there are a variety of states of quantum fields 
which violate them , the most well-known of which is arguably the Casimir vacuum 

I- 

The extent to which quantum field theory allows violations of the local energy 
conditions is not yet completely clear, although progress has been made in recent 
years to answer this question. Two principal approaches for determining the degree 
of violation have been "averaged energy conditions" and "quantum inequalities." The 
first, originally due to Tipler 0, involves a suitable averaging of the local conditions 
over timelike or null geodesies. It can be shown that many of the standard results of 
classical general relativity obtained via global techniques can be proved using only 
the averaged, rather than the local, energy conditions Pj]-[1~0|. Here we take the 
"averaged weak energy condition" (AWEC) to be 

T^u^ifdr > 0, (3) 

where the integral is taken over a timelike geodesic with tangent vector and param- 
eterized by the proper time r. Similarly, we take the "averaged null energy condition" 
(ANEC) to be 

roo 

T M „ K"K v dX>0, (4) 

oo 

where the integral is taken over an affinely parameterized null geodesic with tangent 
vector and affine parameter A There has been a great deal of effort in 



the last several years to determine whether quantum field theory enforces averaged 
energy conditions. Most of this activity has been concentrated primarily on ANEC 



1 12]- [17]. This is in part due to the discovery that violations of ANEC are required 



to maintain traversable wormholes |18, [19). It appears that although ANEC holds 



for a wide class of quantum states in a variety of spacetimes, it does not hold in an 
arbitrary four-dimensional curved spacetime (see Refs. [14], [15], and [16| for more 
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detailed discussions). However, it is quite plausible that a suitable generalization of 
ANEC may hold in more general spacetimes fl?| , [20|| . 

The second approach involves uncertainty principle-type inequalities, derived from 
quantum field theory, which restrict the magnitude and duration of negative energy 
fluxes or densities [f2(], [21], |22| . For example, one such "quantum inequality" (QI) for 
negative energy fluxes seen by inertial observers in two-dimensional flat spacetime 
has the form: 

\F\(Ar) 2 <l, (5) 

where \F\ is the magnitude of the negative energy flux and At is its duration. This 
inequality implies that AE, the amount of negative energy transmitted in time At, 
is limited by 

|A£|At<1. (6) 

Therefore, AE is less than the quantum uncertainty in the energy, (At) -1 , on the 
timescale At. Recently, more precise versions of these inequalities have been derived 
231. These involve an integral of the energy flux multiplied by a "sampling function" , 



i.e., a peaked function of time with a time integral of unity and characteristic width 
t . A convenient choice |2JJ of such a function is t /[7t(t 2 + t 2 )]. If the integrated 
flux, F, is defined by 

then these inequalities may be written as 

1 

F > 

and 



167TT 2 



F on 9 A 1 (9) 

32ii z t 4 

for all To, in two- and four-dimensions, respectively. These inequalities have the 
form required to prevent macroscopic violations of the second law of thermodynamics 
pl| , P3J . It was also discovered that similar inequalities hold for a quantized massless, 
minimally-coupled scalar field propagating on two- and four- dimensional extreme 
Reissner-Nordstrom black hole backgrounds. These inequalities were shown to foil 
attempts to produce an unambiguous violation of cosmic censorship by injecting a 
negative energy flux into an extreme charged black hole |25fl . The latter results 
provide perhaps one of the strongest reasons for the belief that the production of 
large-scale effects via manipulation of negative energy may be forbidden by quantum 
field theory. Classically, any amount of negative energy, no matter how small, injected 
into an extreme black hole is enough to produce a naked singularity. Therefore, if 
one had any chance of producing gross effects with negative energy, the scenario 
discussed above would seem to have offered the best possibility of success. It should 
also be emphasized that the energy-time uncertainty principle was not used as input 
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in the derivation of any of the QI restrictions. They arise directly from quantum field 
theory. 

In a recent paper, we have shown that there exist deep connections between 
averaged energy conditions and quantum inequalities, at least in flat spacetime p3f . 
For a quantized massless, minimally-coupled scalar field in two- and four- dimensional 
Minkowski spacetime, we have derived analogous QI's to Eqs. (§) and (Q) for energy 
density. In particular, for timelike geodesies in two-dimensional flat spacetime, 

to y°° (T„ v u»u")dT 1_ 

7T J-oo T 2 + T 2 ~ 87TT 2 ' 1 J 

for all To- Here (T^u^u") is the renormalized expectation value taken in an arbitrary 
quantum state | ip). In the limit as tq — > oo, we "sample" the entire geodesic, and 
our inequality Eq. (|TCi| ) reduces to AWEC. Note that it is possible that the integral 
in Eq. ([K]) converges, even though the associated AWEC integral in Eq. (|3]) diverges. 
Analogous results hold in this case for null geodesies. 

For a quantized massless scalar field in a two-dimensional Minkowski spacetime 
with compactified spatial dimension, it was also discovered that the difference be- 
tween the expectation values of T^u^u", or T^K^K" , in an arbitrary quantum state 
and in the Casimir vacuum state obey AWEC and ANEC-type inequalities [|20|| . This 
is surprising since it is known (T^j that AWEC and ANEC are violated for (T^u^u") 



and (T^K^K"), respectively, in the (renormalized) Casimir vacuum state by itself. 
Such "difference inequalities" might provide new measures of the degree of energy 
condition violation in cases where the usual averaged energy conditions fail. This ap- 
proach has recently been generalized to arbitrary two-dimensional curved spacetimes 
by Yurtsever [16| . 



For a quantized massless, minimally-coupled scalar field in four-dimensional Minkowski 
spacetime, we also derived the following inequality for timelike geodesies 

to r°° (T^u»u u ) dr 
7i J-oo t 2 + t 2 - 32tt 2 V 



for all To. From this inequality, it was shown that one can derive AWEC and ANEC 
||20|| ■ (The application of QI bounds to constrain the dimensions of traversable worm- 
holes is discussed in a separate manuscript PBJ.) 



In almost all studies of ANEC for quantum fields, the bounds on the ANEC 



integral have been taken from — oo to +00 [12|, [1J, |27j. However, a question of some 
interest is whether ANEC is satisfied in the spacetime of an object collapsing to form 
a black hole. The answer to this question might determine, for example, whether 



Penrose's singularity theorem |28| will still hold in the presence of local violations 



of the energy conditions, such as the Hawking evaporation process p0|| . In this 
case, one is usually concerned with the focusing of null geodesies which generate the 
boundary of the future of a trapped surface. To prove Penrose's theorem, one would 
want ANEC to hold over these half-complete geodesies M, H. This is a stronger 
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condition to impose than to demand that ANEC hold over full (i.e., past and future- 
complete) geodesies. However, even if this version of ANEC fails in some regions 
of a given spacetime, it may hold in enough other regions so that the conclusions 
of Penrose's theorem will still be valid. In addition, it should be noted that ANEC 
is a sufficient condition to insure continued focusing of null geodesies. One can also 
guarantee focusing with other conditions, for example by allowing the right-hand-side 
of the ANEC integral to be only periodically non- negative |7|, or even negative but 
bounded (see Lemma 3 of Ref. §, Eq. (5) of Ref. g], and Ref. Q). 

In the current paper, we examine quantized massless scalar, and electromagnetic 
fields in the Unruh vacuum state in two- and four-dimensional Schwarzschild black 
hole backgrounds. We evaluate AWEC and ANEC integrals along half-complete 
timelike and null geodesies, respectively. Our goal is to determine the regions of 
these spacetimes in which the averaged energy conditions fail, and whether there 
exist any bounds on the extent of the violation. The paper is organized as follows. 
The analysis of 2D evaporating black holes is presented in Sec. (|2|); the results for 
4D black holes are given in Sec. (^). The latter section also includes a discussion of 
ANEC along orbiting null geodesies. Concluding remarks are given in Sec. ([|). Our 
units are taken to be: h = G = c = 1. 



2 2D Evaporating Black Holes 

In this section, we study AWEC and ANEC for a quantized, massless scalar field in the 
Unruh and Boulware vacuum states, along geodesic segments in a 2D Schwarzschild 



spacetime. The discussion here is a natural extension of the analysis in Ref. [31]]. 



2.1 Timelike Observers 

Let be the two- velocity of an inertial observer. The energy density in this observer's 
frame is given by 

U = T^u v . (12) 

In the following discussion, we understand T^ u to denote a quantum expectation value 
in a specified vacuum state. The metric is 

ds 2 = -Cdt 2 + C- 1 dr 2 , (13) 

where C = 1 — 2M/r. An outgoing geodesic observer's two- velocity is 

d4) 

The constant k is the energy per unit rest mass. In our two-dimensional discussion, 
we will consider observers moving in both the Unruh and Boulware vacua. The stress 
tensor components in the Unruh vacuum are B2| : 



_ 1 (7M 2 _4M 1 \ 
Ut ~ 24vA r3 + 32M*J ' {Lb) 
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Ti 



tr 



24vr 



2M\- X 1 
~V) 32M 2 



and 



T 



1 / 2My 2 /M 



24vr V 



32M 2 



The corresponding components in the Boulware vacuum are given by: 



( 



7M 2 AM 



and 



T. 



247TV r 

Ttr = 

1 M 2 
~ 2Aix r 4 



2M 



For an outgoing timelike observer in the Unruh vacuum, 



1 



24vr 



c- 2 k 2 



6M 2 AM 



1QM 2 



+ 



CM 2 



C kVk 2 -C 



32M 2 



16M 2 



(16) 
(17) 

(18) 
(19) 

(20) 
(21) 



Let us consider, for simplicity, the case of an observer who starts out initially very 
close to the horizon and is shot outward at large velocity, i.e., we are interested in 
the limits 

e<M, fc»l, (22) 

where r = 2M + e is the observer's initial position. In this limit k ^> C, so from the 
geodesic equation, Eq. (|I~4D, we have that 



T 



r - 2M - e 

jfe 



+ 0{k- 



(23) 



We now wish to multiply T fJbV u^u l/ by a sampling function, i.e., a peaked function of 
time whose time integral is unity. In Refs. |£0| and [22|, this function was chosen to 
be 

(24) 



7T 



dr 



oo r 2 + r 2 



where the integral was taken over complete geodesies. However, unlike in Ref. [^0] 



here we are integrating over half-infinite geodesies. Since the function given by 
Eq. fl2"4|) is symmetric about r = 0, in our case we may choose 



2r 



dr 



it Jo t z + To 



1 



(25) 



where the proper time r is initialized when the observer starts at r = 2M + e. We 
multiply T fJiU u' 1 u 1 ' by this sampling function, use Eq. (^), and expand the numerator 
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in inverse powers of k. The integral is performed using MACSYMA and the result is 
expanded in powers of e, yielding 



2to 

7T JO 



00 T^uWdr 



k 



To' 



+ 



24tt 2 r 
M 



1 J_ 

e + M n V2M 



+ 0(6° 



In 



kr 



3ir 2 To 3 *"V2M 
Therefore, to leading order in k and e, we have 



0(k 



-is 



2r n A 00 T^uVdr 



'0 

7T JO 



where 



r 2 + r 2 



<5r 



24vr 2 r 5r ' 



e 



(26) 



(27) 



(28) 



Note that in these limits St <C M and 5r <C to 

What is the physical significance of <5r? The negative energy density drops off 
very rapidly with increasing r. In the limit k ^> 1 and e <C M, if we consider the 
proper time spent by the observer in the region of appreciable negative energy to be 



St I r2M+xe \ 

dr ~ — / dr = — t[x 

k J2M+e k 



1) 



(29) 



then typically x — 1 ~ 0(1), so 



5r 



k 



(30) 



To see this a little more explicitly, note that near the horizon, Eq. (|2l|) becomes 

A- 2 



48vr r - 2M 



(31) 



If the observer starts at r = r = 2M + e, then at r = r\ = 2M + 2e, 

fr u^w 1 ^ 1 « - (t 

\- L ^tu Uj " I |r=ri ~ . l^M 



- £tf Lb Lb j \ r =TQ > 



(32) 



so the energy density drops off to 1/4 of its initial value in a distance Ar = e, 
corresponding to a proper time 5t = e/k. Therefore, we should consider the time 
interval 5t not to be the entire proper time over which the energy density is negative, 
but rather a time scale which characterizes a change in the energy density. 
For an outgoing timelike observer in the Boulware vacuum: 



» vU U 24vr 



6M 2 4M 



+ 



CM 2 



(33) 
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Note that this quantity is negative everywhere for r > 2M. If we perform the same 
set of operations as discussed above, we find that Eq. (p7|) holds for the Boulware 
vacuum as well. In the r — > oo limit of Eq. fl2T|), we "sample" the entire (i.e., 
half-infinite) geodesic and find 

T^uWt « --i- , (34) 
o Aottot 

and hence 

roo | 



. T^u^dr >-—■ (35) 
/o or 

The dominant contribution to the integral in Eq. (0) is from the lower limit, since 
Tf JLV u^u v in the Boulware vacuum state also has the form given in Eq. (|3lD near the 
horizon, but drops off rapidly with increasing r. By contrast, the tq — > oo limit of 
Eq. (p7[) for the Unruh vacuum state results in an integral which diverges positively 
at large r. This is due to the fact that in 2D the Hawking radiation does not drop 
off with distance, but remains constant at large r. Therefore, AWEC is satisfied in 
the Unruh vacuum state in 2D. However, as discussed in Sec. (0), the situation is 
different in 4D. We could perform the integral on the left-hand- side of Eq. (|3~4| ) for 
the Unruh vacuum state, but instead truncate the integration at some finite value of 
r = r max , with 2M + e < r max < oo. For fixed r max , by making e small enough we 
can always arrange it so that the negative energy contribution from the lower limit, 
given by Eq. (|3~I|), dominates the integral. In that case, using Eq. ([51]) , we have that 

T^uVdr^-——>--. (36) 
o 487TOT or 



The inequalities, Eqs. ( |35| ) and fl36|), represent bounds on the degree of AWEC 
violation seen by timelike geodesic observers who start out very close to the horizon 
and are shot outward at high speed. In this limit, the longer the (proper) timescale 
over which the observer sees a significant change in the energy density, the smaller 
is the magnitude of the integrated negative energy density seen by that observer. 
Since the negative energy density drops off rapidly with increasing r, to remain in 
the negative energy density the observer must stay close to the horizon. The closer 
the observer is to the horizon, the larger is the magnitude of the negative energy 
density. However, in order to remain close to the horizon for a long time as seen 
by a distant observer, the observer's trajectory must be nearly lightlike. Therefore, 
although the observer spends a long time in the negative energy region as seen by the 
distant observer, the proper time spent in the region of appreciable negative energy 
(as measured by St) decreases with the observer's proximity to the horizon. These 
inequalities are similar in form to Eq. @. 



9 



2.2 Null Geodesies in 2D: Off the Horizon 



We now wish to examine ANEC for null geodesies in 2D Schwarzschild spacetime. 
From Eq. flTB|), and from the equation for null geodesies, we have 



— = ±C- l E, — 
dX dX 



±E. 



(37) 



where E is an arbitrary positive constant, whose value fixes the scale of the afline 
parameter, A. We want to examine ANEC along both ingoing and outgoing future- 
directed (i.e., dX > 0) null geodesies, inside and outside the horizon. There are four 
cases (see Fig. 1): 



E > , for r > 2M (outgoing rays) , (38) 

—E < , for r > 2M (ingoing rays) , (39) 

—E < , for r < 2M (outgoing rays) , (40) 

—E < , for r < 2M (ingoing rays) . (41) 



for future-directed null rays. 

If we use Eqs. ([38|) - (|4l|), and Eqs. (^) - (|17D, we obtain for the Unruh vacuum 
state 



1) 


dt 
dX 


> 


o, 


dr 

dX = 


2) 


dt 
dX 


> 


o, 


dr 

dX = 


3) 


dt 
dX 


> 


o, 


dr 

dX = 


4) 


dt 
dX 


< 


o, 


dr 

dX = 


Note that C = 


= (1 


-2M/r) 


chan 



E 2 

^ 24vr 



1 



2M 



-2 



6M 2 AM 



and 



E^ 

247T 



2M 



6M 2 



4M 



+ 



(outgoing) , (42) 



(ingoing) . (43) 



8M 2 



To obtain Eq. fl43|), we have used the fact that the sign of T tr changes inside the 
horizon, as can be seen from Eq. (|16|), together with Eqs. (|39D and (|4lD . Note that 
for Eq. (|43D, the term in square brackets vanishes as r — > 2M. A Taylor expansion 
shows that, to leading order, it vanishes as (r — 2M) 2 , so it will cancel the (r — 2M) 2 
divergence in the denominator. Therefore Eq. (f43D is finite on the horizon, as it should 
be. However, there appears to be a discontinuity in Eq. ( ^2| ) for an outgoing null ray 
on the horizon, r = 2M. The situation for null rays on the horizon is sufficiently 
subtle as to warrant a separate discussion. This is provided in the next subsection. 

From Eq. ([4*2]), T^ lu K fl K l/ is negative for outgoing null geodesies when 
r > 1.5M. Therefore the local null energy condition, Eq. (Q) is violated along outgo- 
ing null vectors slightly inside the horizon, as well as along all outgoing rays outside 
(and on) the horizon. For r > 2M, consider an outgoing null geodesic starting at 
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r = 2M + e, with e <C M. Examine ANEC along this ray to obtain 



h = 



T llu K f "K u dX 



E 

24:71 J2M+ 



(r - 2M 



6M 2 4M 



(44) 



where we have used Eq. fl38|). Here A = is the value of the affine parameter at 
r = 2M + e. Equations (^) and (j44[) apply for the Boulware vacuum as well as for 
the Unruh vacuum. This is due to the fact that in the latter case, the outgoing null 
geodesies are parallel to the Hawking radiation. Hence the terms which would cause 
the integral in Eq. Q3BD to diverge as r max — > 00 are absent in the null geodesic case. 

For r < 2M, along an outgoing null geodesic starting just inside the horizon at 
r = 2M-e: 



h = 



T^K^^dX 



E 

2M 



2M-e 



2M 



6M 2 4M 



dr . 



(45) 



where we have used Eq. (fiOD , and A/ is the value of the affine parameter at r = r min , 
the minimum value of r attained by the null geodesies. Notice that the integrand in 
Eq. (|45|) will be dominated by the positive 6M 2 /r 2 term which diverges for small r. 
From Eq. (|43D , we see that T^ U K^K U is positive everywhere for ingoing null geodesies. 
Thus we immediately see that ANEC is satisfied for these two sets of null geodesies. 

In Eq. we perform the integration and then expand the result in the limit of 
small e to find 

hm h ~ + 0(log(e/M)) < . (46) 

e^O 4S7T e 

Note that 5X = E/e is the characteristic affine parameter distance over which the 
negative energy is decreasing along the outgoing geodesic. We can rewrite Eq. (f46|) 

as 

roo 1 1 







T^K v dX 



48tt(5A 

Note that it is invariant under rescaling of the 



This is the null version of Eq. 
affine parameter. 

Similarly, for Eq. (f4"o]), we may perform the integration and expand the result in 



(4£ 



the small r m ,„ limit. For fixed e, the result is 



E 

lim I 2 ~ 



+ 0(ln{r min /M)) > 0. 



However, for fixed r min and e — > 0, -Z2 has the same behavior as I± : 



487T6 



+ 0(ln(e/M)) < 0. 



(49) 



Thus we see that ANEC is violated for outgoing null rays just outside the horizon, 
but it is satisfied for outgoing null rays inside the horizon in the limit r min — > 0. (The 
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divergence in Eq. ( pE6| ) as e — > may be circumvented by an appropriate choice of 
scaling of the affine parameter, i.e., by an appropriate choice of E. This and related 
issues will be discussed in the next sub-section.) Outgoing rays which originate in 
the region 1.5M < r < 2M initially encounter negative energy, and later encounter 
positive energy near the singularity, r = 0. This is analogous to the case of outgoing 
timelike geodesies in the Unruh vacuum, which encounter negative energy near the 
horizon and positive energy at large distances due to the Hawking radiation. We also 
see that ANEC holds for ingoing null geodesies, whether they originate outside or 
inside the horizon. 

One might at first sight conclude that the reasoning here is a bit circular. It 
could be argued that it is no surprise that evaluation of the ANEC integral for null 
geodesies inside the horizon yields a divergent positive result. After all, we placed 
our quantized field on a background which already had a singularity in it. Thus one 
could argue that divergent values of our ANEC integrals for these geodesies are not 
unexpected, because the behavior of the fields inside the horizon is dominated by 
the singularity. However, it is not obvious, a priori, what the sign of these integrals 
should be. Since these are test fields on a given background, as opposed to fields 
which generate the background spacetime (which we cannot have in any case in 2D), 
it is possible that the integrals might have turned out to diverge with either sign. 
Similar ANEC integrals for quantized massless fermion fields in 2D have the same 
signs as for the massless scalar field. This follows from the fact that the renormalized 
fermion stress tensor is identical to that in the scalar case Pl |, even though the 
formally divergent tensors have opposite signs in the two cases. 

2.3 Null Geodesies in 2D: On the Horizon 

Recall that our results for outgoing null geodesies outside the horizon are formulated 
in terms of the null tangent vector 

K» = E(C-\ 1). (50) 

where the components are given in terms of Schwarzschild (t, r) coordinates. Since 
these coordinates are badly-behaved on the horizon, let us switch to Kruskal null 
coordinates, given by 

U = -e~ KU , V = e KV , (51) 

where u = t — r*,v = t + r*, with r* the usual tortoise coordinate and k = 1/(4M). 
In these coordinates, our null vector has the form 

M§.£)-( "''")■ 

If E has the same value for all outgoing null geodesies, then would not be defined 
in the limit where the geodesies approach the horizon, since C — > 0. However, E need 
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only be constant along each null geodesic, but not the same constant for different null 
geodesies. 

Consider a sequence of outgoing null geodesies which start at different values of 
r, e.g., the histories of a sequence of photons emitted by an infalling observer. Let tq 
be the value of r at which the null geodesic begins (which will be different for each 
null geodesic). We might choose the various values of E for different null geodesies 
by the following prescription. Let the observer emit photons of fixed frequency u e in 
his rest frame. The time-component of the observer's two-velocity (in Schwarzschild 
coordinates) is 

, dt ( 2MY 1 - , , 

" = 3; = I 1 — J E < < 53 > 

where E is the observer's energy per unit rest mass, measured at infinity. This is 
just the rate of a clock at infinity as compared to the rate of the clock carried by the 
infalling observer. The frequency of a photon at infinity is u^, which is related to u e 
by 

«. = dr = t _ 2JfU, 
uj e dt \ r J 

Note that the expression on the right-hand-side includes both the Doppler effect and 
the gravitational redshift. The factor of E~ l reflects the fact that the faster the 
observer is shot toward the black hole (i.e., the larger E is), the greater will be the 
Doppler shift. If we choose a scaling for the affine parameter such that E = u^, for 
fixed u e , then 

E = (l - u e E~ l oc C(r ) . (55) 

Now E/C{tq) is finite in the limit r$ — > 2M, so that the components of are finite 
on the horizon, in either Schwarzschild or Kruskal coordinates. 

However, is still not the same as k^, the affinely parameterized null tangent 
vector on the horizon. In Kruskal null coordinates the vector has components 

k p = (0, k) . (56) 

(See the discussion on p. 331 of Ref. |]35||.) On the horizon, 




K , = er(0, K ), (57) 



=2I\1 



so these vectors still differ by a factor of exp(nv). Thus, a null vector such as K^, 
which is an affinely parameterized null tangent vector off the horizon, does not nec- 
essarily remain affinely parameterized on the horizon. 

With these lessons in hand, let us now consider the integral / T^p kP 1 k v dV along 
a portion of the future event horizon. Here Tp_ p is the vacuum expectation value of 
the stress-tensor in the Unruh vacuum state expressed in Kruskal null coordinates, 
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and kP 1 is an affinely parameterized null tangent vector on the horizon. The Kruskal 
advanced time coordinate V, is an affine parameter on the horizon (see for example, 
p. 331 of Ref. |35| or p. 122 of Ref. ||36|| ). From Eq. fl5"6|), and a straightforward 



coordinate transformation of Eqs. (p~5|)- (|ITD , we obtain 



T QP k 11 k v dV = — . (58) 

v 768irM 2 V 1 ; 

This integral is taken over a portion of the future horizon; Vq is the value of the 
V = const line which intersects the future horizon at the (arbitrary) event where we 
start our null geodesic. 

The right-hand-side of Eq. ( |5*5D goes to — oo as Vq — > 0. However, this would 
correspond to integrating along the full future horizon of an eternal black hole space- 
time. The Unruh vacuum state is singular on the past horizon (i.e., V — 0) of such 
a spacetime, so it is not surprising that we get a divergence in this case. Physically, 
this would correspond to a black hole that Hawking-radiates for an infinite time. For 
the Unruh vacuum state to be realizable, we must consider only physically realistic 
collapse spacetimes in which the black hole forms at a finite time in the past. There- 
fore, the smallest allowed value of V = V min , which is the value of the V = const 
line at which the the worldline of the surface of the collapsing star intersects the 
future event horizon, as shown in Fig. (2). For such spacetimes, the ANEC integral 
along a portion of the future horizon is negative but finite. This result seems to be 
in agreement with the results of Wald and Yurtsever |14|. They find that ANEC is 
satisfied along the complete null geodesic comprised of the future horizon plus the 
null line which would have been the past horizon in an eternal black hole spacetime. 
In a collapse spacetime, this complete null geodesic originates at past null infinity, 
propagates through the collapsing star, "reflects off" r = 0, and then emerges into 
the vacuum outside the star at Vo = V m i n (see Fig. (2) ). Evidently there is a pos- 
itive contribution to the ANEC integral from the part of the geodesic prior to its 
exit from the collapsing star ||37[ . This contribution must be greater than or equal 
to (768 7r M 2 Vmin)' 1 . (Note that the source of the positive contribution we refer 
to here is the quantized field, not the classical collapsing matter. The latter would 
presumably make the ANEC integral even more positive.) 

The dimensions of kP are (length) 1 , whereas Tpp, V, and Vo are dimensionless. 
Let us now rescale the coordinates to get a more familiar choice of dimensions. Let 
x d = i-e . ; (j = K -iU,V = k~ 1 V, so V = k~ x Vq. Then = n 2 T pp , and 

k^ = k' 1 kP 1 . If we rewrite Eq. ([58]) in terms of the rescaled coordinates, we obtain 

T r ,o P kP dV ' 



Vo " V 7Q8ttk 2 M 2 V 



o 

' (59) 



48tt V 



o 



Compare this result with our result for null rays outside the horizon, i.e., Eq. fl4"7j). 
in which 5X is interpreted as the characteristic affine parameter distance along the 
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null geodesic over which the negative energy density falls off. Let us now show that 
a similar interpretation holds on the horizon. The integrand of Eq. (p9|) is equal to 
-1/(48 7i V 2 ), so 



oo „„„ 1 foo I I \ 

T^WWdV = / \^-)dV 

1 



487Tl/n 



(60) 



Note that V$ can be interpreted as the characteristic affine parameter distance over 
which the integrand falls off, i.e., as V increases from V = Vo to V — 2Vq, k v 
falls off to 1/4 of its initial value. Thus the interpretation of the QI for null geodesies 



on the horizon is the same as that for Eq. (47) 



3 4D Evaporating Black Holes 

In this section, we examine AWEC and ANEC in four-dimensional evaporating black 
hole spacetimes. In these spacetimes, the components of the renormalized expectation 
value of the stress-energy tensor are only known numerically. Accordingly, we make 
use of the numerical data of Elster |38) for the conformally-coupled massless scalar 
field and that of Jensen, McLaughlin, and Ottewill |3!| for the electromagnetic field, 
evaluated in the Unruh vacuum state. We examine the cases of outgoing radial 
timelike observers, outgoing radial null geodesies, and orbiting null geodesies. Only 
geodesies outside the horizon are considered, as no numerical data are available for 
r < 2M. We find that although AWEC and ANEC are violated in all three cases, 
there appear to be Ql-type bounds in each case which constrain the degree of the 
violations, as in two-dimensional spacetime. 

These results might at first sight seem to be in disagreement with recent results 
of Visser E0[ . He found that for the renormalized vacuum expectation value of the 



stress-energy tensor of a conformally-coupled massless scalar field in Schwarzschild 
spacetime, ANEC is satisfied for all null geodesies that reach spatial infinity, and for 
orbiting geodesies at r = 3M . However, his results are applicable for the Hartle- 
Hawking vacuum state, which represents a black hole in thermal equilibrium with a 
surrounding external thermal bath of radiation. It is therefore not too surprising to 
find that ANEC holds in that case. By contrast, our analysis is performed in the 
Unruh vacuum state, which corresponds to the situation of a black hole evaporating 
into empty space. 

3.1 Outgoing Radial Timelike Observers 

The four-velocity of an outgoing radial timelike observer in 4D Schwarzschild space- 
time is given by 

u" = (A Vf^c, 0, 0) , (61) 
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where again C = (1 — 2M/r). Let us adopt the notation of Elster [RS so that: 
/i(r) = — T l t , p r (r) = T r r , L = —Attt 2 C T tr , where the components of the stress- 
tensor here represent renormalized expectation values in the Unruh vacuum state. 
Therefore, we may write 

k 2 , kL 



T^vTu" = - (jl + Pr ) - Pr - ^ r2c2 ^[W~C . (62) 

If a Ql-type bound of the kind found in 2D, i.e., Eq. (|35|), also holds in 4D, then we 



might expect it [|4|] to have the form 



oo ] 

r^u'drjfc-— , (63) 

where again St is the characteristic proper time which the observer spends in the 
negative energy. (Note that since the fields we will be considering are massless, the 
mass M of the black hole is only natural length scale which appears in the problem.) 
We want to numerically compare the left and right-hand-sides of Eq. (|63|) , in order 
to see if the inequality is satisfied. 

The integrand in Eq. ( |63"D drops off rapidly with increasing r near the horizon 
(unlike in the 2D case), hence the dominant contribution to the AWEC integral 
comes from the initial part of the trajectory. We will therefore approximate the 
AWEC integral as 

POO 

/ T^uWdra T(r)\ ro ST, (64) 
Jo 

where T(r) = T tiv u^u v is understood to be evaluated at the starting point r (min- 
imum r) of the trajectory. Since the magnitude of T(r) rapidly decreases near the 
horizon, our estimate of using T(r ) gives a result which is somewhat more negative 
than the true result. Therefore, correction of this estimate only goes in the direction 
of strengthening the inequality. 

Our results for the electromagnetic field are displayed in Table 1. The values of 
p, and p r for the electromagnetic field are from the numerical data given in Ref . [^UJ . 
and the value of L = 3.4 x 10~ 4 /M 2 used in this case is the photon luminosity given 
by Page [|42|| . In the fourth column of Table 1, Sr is chosen such that the magnitude 
of T(r) decreases by at least a factor of two from r to r + Sr. The characteristic 
proper time, St, is computed from 

ST = L (65) 

The next to last column in the table is our approximation of the AWEC integral given 
in Eq. fl6"4|), while the last column represents the right-hand-side of our QI given by 
Eq. (|63|). By examining the last two columns of Table 1, we see that the QI, Eq. (|63|), 
is easily satisfied. A similar analysis for the conformally-coupled massless scalar field 
shows that the QI is easily satisfied in that case as well. Numerical values of n,p r , 
and L = 7.44 x 1Q~ 5 /M 2 for the scalar field were obtained from the numerical data 
given in Ref. 
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-5 










10 


0.10 
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1 
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-0.45 
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0.40 


-6.8 x 10~ 3 


-2.5 










10 


0.20 


-1.3 x 10" 2 


-5 



Table 1: Outgoing timelike geodesies moving through the quantized electromagnetic 
field. 



3.2 Radial Null Rays 

The components of the tangent vector to an affinely parameterized outgoing radial 
null geodesic (outside the horizon) are given by 

K» = E(C- 1 , 1,0,0). (66) 

Thus we have 

T^JC = - E ' 2 Unr (r - 2M) [fi(r)+p r (r)\ - l) . (67) 
2tt (r — 2M) { ) 

The 4D analog of the QI for null geodesies in 2D given by Eq. ([47D is 



1 

T^K"d\>-j^, (68) 

where 5X is again interpreted as the characteristic affine parameter distance along 
the null geodesic over which the negative energy density falls off. If we define 

T = {T^K") , (69) 

and use the fact that dr/dX = E and 5 A = 5r/E, then it is easily seen that the QI, 
Eq. fl68D is equivalent to 

£ Tdr *-WT r - < 70) 

where r > 2M is the initial starting value of r for the outgoing null geodesic and Sr is 
the characteristic interval in r over which the negative energy falls off. We numerically 
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Mr )M 4 


p r (r )M 4 


5r/M 


T5r M 3 


-M/8r 
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-3.4 x 10" 5 
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-3 x 10- 7 
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-4.8 x 10- 6 


-0.3 


5 


-1.5 x 10" 6 


-3 x 1(T 7 
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-7.8 x 10" 6 


-0.2 



Table 2: Data is given for outgoing null geodesies, where /j, and p r are the energy 
density and pressure for a massless scalar field. 



r /M 


/i(r )M 4 


p r (r )M 4 


5r/M 


T5r M 3 


-M/Sr 


2.1 
2.2 
2.5 

3 


-2.7 x 10" 3 
-1.8 x 10~ 3 
-1.1 x 10~ 3 
-5.4 x 10" 4 


1.9 x 10- 3 
1.5 x 10~ 3 
6.8 x 10" 4 
3.2 x 10" 4 


0.5 
1 
1 
2 


-0.01 
-0.004 
-0.002 
-1.4 x 10" 3 


-2 
-1 
-1 
-0.5 



Table 3: Data is given for outgoing null geodesies, where /x and p r are the energy 
density and pressure for the electromagnetic field. 



evaluate each side of Eq. ( ffO|) to see if the inequality is satisfied. Following a similar 
procedure to that used for outgoing radial timelike observers, we approximate the 
ANEC integral by 

roo 

/ T dr « T(r ) 6r . (71) 

J To 

As discussed in the last section, this approximation gives a result which is somewhat 
more negative than the true result. Therefore, if the QI holds for our approximation, 
then the inequality in the actual case is even stronger. 

Our results for the scalar field are given in Table 2, and those for the electromag- 
netic field appear in Table 3, where as before the numerical data of Refs. |J8] and 



], were used for the scalar and electromagnetic fields, respectively. The next to last 
column in each table is our approximation of the ANEC integral given in Eq. (|7T|). 
while the last column represents the right-hand-side of our QI given by Eq. (|70|). By 
examining the last two columns of each table, we see that the QI, Eq. (|70|), is easily 
satisfied and hence so is Eq. (|6"8D . Note that the inequalities in 4D seem to be satisfied 
by a wider margin than the corresponding QI's in 2D. This is due to the fact that the 
numerical factors which appear in the components of the renormalized expectation 
values of the stress-tensor are characteristically smaller in 4D than in 2D spacetime, 
at least for radial null geodesies in the scalar field case. 

In our calculations, we have ignored the backreaction of the Hawking radiation 
on the background spacetime. It might be better to repeat our calculations on an 
evolving background, such as a Vaidya spacetime. However, we can estimate the 
effects of backreaction on an outgoing null ray just outside the horizon in the following 
way. Our approximation should hold as long as the time of escape of a null ray from 
near the horizon is small compared to the evaporation time of the black hole. Here 
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these times are measured by an observer at infinity. For such a null ray we have 



cxM 
l2M+e 

Solving for t esc yields 



/ dt = (f -2M/rY l dr. (72) 

J J2M+t 



t esc = xM + 2M\n(^ -lj - 2M - e + 2Mln^^ , (73) 

where t esc is the time for the ray to go from r = 2M + e to r = xM. For i>2 and 
finite, in the limit of small e: 

t esc ^2M\n(^pj . (74) 

By contrast, the evaporation time of the hole is t evap = AM 3 , where A 10 4 /n, 
and n is the number of species of particles, each of which is assumed to contribute 
approximately 10~ 4 M~ 2 to the black hole's luminosity Our approximation should 
break down only when e < 2M exp[-AM 2 /2]. 



3.3 Orbiting Null Rays 

We now examine orbiting null geodesies in Schwarzschild spacetime. The line element 
is 

ds 2 = -Cdt 2 + C~ x dr 2 + r 2 [dQ 2 + sin 2 ^0 2 ) , (75) 

where again C — (1 — 2M/r). For an orbiting null ray in Schwarzschild with 6 = 

n / 2 = const, 

dt_ _r_ d4_ 

d\ d\' { ' 

Here we choose d<p/dX > 0, for our future-directed ray. From the geodesic equations, 
we have 

1 s (77) 
= Tx = 7 2 ' (78) 

where E and L are constants. It follows that 

^ = rC~W. (79) 



At r = 3M, we have 

L 
E 



= 3V3M. (80) 

r=3M 
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From the numerical values in Figs. (1) and (3) of Ref. []39] for the quantized electro- 
magnetic field, and Eq. (ISO), we find that 



7ViW| r=3 M = 3E 2 



-Tt + t; 



\r=3M 



E 2 

w -3.1 xlO" 4 — , (81) 

where = dx^/dX is the tangent to the null geodesic. For one orbit, A0 = 2n = 
(L/9M 2 ) AA, so 

m""aT- (82) 

We may therefore rewrite Eq. (|8T|) as 

T, u K"K u \ r - 3M ss -3.7 x 10~ 2 ? . (83) 

1 M 2 (AA) V ; 



Integration of T^ U K^K V over one orbit gives 



AA 1 

„ T r JC'ICdX«^.1xVf wm . (84) 



Thus we see that for orbiting null geodesies 

Jo """" " M 2 (AA) 



T^K v dX > ~ 775~7a~T\ • (85) 



The expressions given by Eqs. (|8J) and ( p^) are of the same form as those 
for radial null geodesies discussed earlier. The difference here is that we integrate 
only over one orbit, i.e., we do not consider multiple traversals through the same 
negative energy region. The inclusion of the contributions of an arbitrary number 
of orbits would violate the inequality given in Eq. (85|). (Since the r = 3M orbit 



is unstable, realistically we would expect the null ray to at best only make a few 
orbits.) However, it is of interest to note that if we consider the null geodesies which 
generate the boundary of the future of an event P with r = 3M, then the orbiting 
null geodesies leave this boundary after half an orbit. This is due to the fact that 
orbiting null geodesies originating at P, which initially set off in opposite spatial 
directions, encounter a conjugate point (i.e., the geodesies cross) on the far side of 
the black hole after half an orbit P3]| , f4*4}| . At this point, the geodesies leave the 



boundary of the future of P (i.e., an event Q which lies beyond the conjugate point 
along either null geodesic can be connected to P by a timelike curve), and therefore 
are not achronal. 

This observation may be relevant in the context of proving singularity theorems 
in the presence of negative energy. For example, in the proof of a singularity theorem 
such as the Penrose theorem, one uses either the local or the averaged null energy 
condition to insure the focusing of the null geodesies which generate the boundary of 
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the future of a closed trapped surface. Geodesies of the orbiting variety considered 
here would not relevant, at least in the spherically symmetric case, since they leave 
the boundary of the future of a point in any case, irrespective of energy condition- 
focusing arguments. We also note in passing that, at least in this simple case, while 
these null geodesies are in the boundary they obey the inequality, Eq. ( |35*D . 

Let us compare the previous results with the following situation. Consider a 
4D Minkowski spacetime with one spatial dimension compactified, and a quantized, 
minimally-coupled, massless scalar field in the Casimir vacuum state on this back- 
ground. Let K 1 * = E(l, 1,0,0) be the tangent to a null geodesic in the x-direction. 
The vacuum expectation value of the stress-tensor in this state is 

T liV = -^j xdiag {-1,-3, 1,1}, (86) 

where we have taken x as the compactified dimension with circumference L. It follows 
that 

4ti 2 E 2 

T^K»=-—. (87) 

The affine parameter lapse for one orbit is X p = L/E. If we integrate over one orbit, 
we find that 



f Xp T^KWdX = (t^K»K»)(\ 

J 



4vr 2 



45 L 2 X p 



We see that this result has the same form as Eq. ( p4|) for the black hole case. Both of 
these expressions are invariant under rescaling of the affine parameter (as is Eq. (|6S|) 
for radial null geodesies). Here, as in the black hole case, we have ignored the possi- 
bility of multiple orbits, which would violate the inequality. However, it is possible 
that in the 4D compactified case, as was found previously in 2D ||20| , an ANEC-type 



bound might hold on the difference between the expectation values of T flu K^K u in 
an arbitrary quantum state and in the Casimir vacuum state. 



4 Conclusions 

In this paper, we have examined AWEC and ANEC in both two and four- dimensional 
evaporating Schwarzschild black hole spacetimes, for quantized massless scalar and 
electromagnetic fields. Our 2D results are the following: It was discovered that 
AWEC is satisfied for outgoing radial timelike geodesies which start just outside the 
horizon and reach infinity, for the Unruh vacuum state. This is due to the positive 
contribution to the AWEC integral of the Hawking radiation at large distances, which 
outweighs the contribution of the negative energy near the horizon. By contrast, 
AWEC is violated for such geodesies in the Boulware vacuum state. Similarly, ANEC 



21 



is violated for half-complete outgoing radial null geodesies, which start just outside 
the horizon, in both the Unruh and Boulware vacuum states. This is also true on the 
horizon in the Unruh vacuum state. However, ANEC is satisfied along ingoing null 
geodesies and along outgoing null geodesies inside the horizon in the Unruh vacuum 
state, due to a large positive contribution to the integral near the singularity. 

In 4D, we examined AWEC and ANEC in the Unruh vacuum state, using pre- 
viously calculated numerical values of the stress-tensor components. There it was 
found that AWEC is violated for outgoing timelike geodesies which start just outside 
the horizon and reach infinity. (Note that in 4D, the Hawking radiation decays like 
1/r 2 , so the large positive contribution to the AWEC integral found in 2D does not 
arise here.) Similarly, we found that ANEC is violated for analogous outgoing null 
geodesies. This condition is violated for orbiting null geodesies as well. 

However, in all cases where the conditions are violated, there appear to be Ql-type 
bounds which limit the degree of the violations. Let St be the characteristic proper 
time for timelike geodesies, and SX be the affine parameter length for null geodesies, 
over which the negative energy density changes significantly. In 2D the QI bounds 
have the form that the integrated negative energy density along the timelike or null 
geodesic (i.e., the AWEC or ANEC integral) is greater than or equal to minus the 
inverse of St or SX, respectively. In 4D, the corresponding inequalities have additional 
factors of M ~ 2 on the right-hand-sides. 

Can the semi-classical effects of negative energy, in processes such as black hole 
evaporation, invalidate the singularity theorems before quantum gravity effects be- 
come significant? Although we do not have a definite proof in the general case, our 
results suggest that the answer may be no. It remains possible that even though 
ANEC (or even a weaker energy condition [17]) might fail in some regions of a given 



spacetime, it may hold in enough other regions for the conclusions of the singularity 
theorems to be maintained. Our (admittedly) 2D results for ANEC would seem to 
support this view. In Ref. [[J, it was shown that if ANEC is satisfied along the null 
geodesies which generate the boundary of the future of a trapped surface, then Pen- 
rose's singularity theorem will still hold. There the local null energy condition was 
replaced by ANEC along half-complete null geodesies. However, two points should 
be emphasized about that result. First, it is not necessary to assume that the ANEC 
integral is non- negative to prove the result 0, |l7j. In fact, from Lemma 3 of Ref. || 
or Eq. (5) of Ref. |J and the Einstein equations, it is sufficient to have 

T, v K»K»dX> ^, (89) 

^ 

where 9 < is the initial expansion of the null geodesic congruence at the trapped 
surface. The ANEC integral may be negative as long as it is not more negative 
than 6q/2. Physically this implies that a singularity will still form provided that 
the defocusing effects due to the presence of any negative energy are more than 
compensated for by the initial convergence of the null geodesies produced by the 
trapped surface. 
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Second, the existence of only one such trapped surface whose orthogonal null 
geodesies obey either ANEC or the weaker inequality, Eq. fl89|), is required to prove the 
occurrence of a singularity. There may exist other trapped surfaces in the spacetime 
whose orthogonal null geodesies do not obey either of the inequalities. However, 
as long as there exists at least one trapped surface with the desired properties, a 
singularity is inevitable. 

Consider an evaporating black hole spacetime, with backreaction included. The 
Penrose diagram for the standard scenario is shown in Fig. (3). Let us assume that 
this standard picture is correct. The line H + is the event horizon, the dashed curve 
A is the apparent horizon, and the lines gi, g 2 , and g 3 are outgoing null geodesies 
orthogonal to the trapped surfaces labelled 1, 2, and 3, respectively. The outgoing 
null rays from 1 do not focus but instead reach future null infinity, so one expects 
that any averaged energy condition which would insure focusing, such as ANEC or 
Eq. (|39D , would fail along them. The closest analogs to these rays in our analysis 
are the outgoing null rays just outside the horizon shown in Fig. (1). It was found 
that in both two and four- dimensional Schwarzschild spacetime, ANEC was violated 
along these rays, although in each case a QI limits the degree of the violation. One 
would expect that if an analogous QI holds in more general cases, such as Fig. (3), 
it will not be strong enough to guarantee focusing of null rays in the region between 
H + and A in Fig. (3). We could not examine the weaker inequality Eq. (|89|), since 
in a static evaporating black hole background there are no trapped surfaces outside 
the horizon. For the outgoing null rays from 2 in Fig. (3) (a trapped surface inside 
the collapsing matter), one would suspect that ANEC probably holds. In this case, 
the focusing effects of the matter of the star most likely overwhelm any defocusing by 
negative energy. For the outgoing null rays from 3 (a trapped surface "behind" the 
event horizon, in the vacuum outside the collapsing matter), one might expect either 
ANEC or a similar but weaker inequality such as Eq. ( |5§D to hold, since the null rays 
eventually do run into a singularity. This is assuming, of course, that the "standard" 
scenario depicted in Fig. (3) is in fact the correct one - a view which could well be 
wrong. In our 2D analysis, we found that ANEC was satisfied along these rays, and 
also along ingoing null rays. 

We showed that our approximation of neglecting backreaction should fail only 
when the time of escape of a null ray from near the horizon is comparable to the 
evaporation time of the black hole. This implied that such a ray would have to 
start out extremely close to the horizon. Nevertheless, it would be useful to attempt 
a similar analysis to ours for evaporating black hole spacetimes with backreaction 
taken into account. The situation is complicated by our ignorance of the detailed 
form of the stress-tensor components in such cases. 

Recent work of Kuo and Ford |45|, [46| has shown that negative energy densities 
in flat spacetime are subject to large fluctuations. In these circumstances, one does 
not expect the semi-classical Einstein equations to be a good approximation. If the 
same is true of negative energy densities in curved spacetime, and we see no reason to 
believe otherwise, then the question arises as to whether such fluctuations could cause 
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a gross failure of the semi-classical approximation at scales well above the Planck 
scale. Large stress-tensor fluctuations would presumably produce large fluctuations 
in spacetime geometry. Therefore, it is very important to determine whether such 
fluctuations might render suspect the semi-classical picture of the Hawking process, 
or the conclusions of the singularity theorems in the presence of negative energy. On 
the other hand it could turn out, for example in the Hawking process, that if the 
timescale of the stress-tensor fluctuations is very short compared to other relevant 
time scales, such as the evaporation time of the hole, then the standard analysis may 
still be valid. In the proofs of the singularity theorems, one assumes a fixed classical 
spacetime background which determines the causal structure. If the background 
geometry fluctuates, it is not clear that the usual global techniques employed in 



classical general relativity will still be useful ||47 |. These questions are currently 
under active investigation. 
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Figure Captions 



• [1] Radial null geodesies in the spacetime of a black hole formed by gravitational 
collapse: an ingoing geodesic (a), an outgoing geodesic inside the horizon (b), 
and an outgoing geodesic outside of the horizon (c). 

• [2] The spacetime of a black hole formed by gravitational collapse. The null 
ray which lies in the future horizon enters the collapsing body (shaded region), 
and then re-emerges into the vacuum spacetime at V m i n . 

• [3] The spacetime of an evaporating black hole, including the effects of backre- 
action. 
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